A POLYHEDRAL APPROACH TO COMPUTING BORDER BASES 



GABOR BRAUN AND SEBASTIAN POKUTTA 

Abstract. Border bases can be considered to be the natural extension of Grobner bases that have several 
advantages. Unfortunately, to date the classical border basis algorithm relies on (degree-compatible) term 
orderings and implicitly on reduced Grobner bases. We adapt the classical border basis algorithm to 
allow for calculating border bases for arbitrary degree-compatible order ideals, which is independent from 
term orderings. Moreover, the algorithm also supports calculating degree-compatible order ideals with 
preference on contained elements, even though finding a preferred order ideal is NP-hard. Effectively we 
retain degree-compatibility only to successively extend our computation degree-by-degree. 

The adaptation is based on our polyhedral characterization: order ideals that support a border basis 
correspond one-to-one to integral points of the order ideal polytope. This establishes a crucial connection 
between the ideal and the combinatorial structure of the associated factor spaces. 



1. Introduction 

Grobner bases are fundamental tools in commutative algebra to model and perform important oper- 
ations on ideals such as intersection, membership test, elimination, projection, and many more. More 
precisely, the Grobner bases framework makes these operations computationally accessible allowing to 
perform actual computations on ideals. 

1.1. Comparing border bases and Grobner bases. Unfortunately, Grobner bases are not always 
well suited to perform theses operations, in particular when the actual ideals under consideration are 
inferred from measured data. In fact, Grobner bases do not react smoothly to small error in the input 
data. Border bases are a natural generalization of Grobner bases that are believed to deform smoothly 
in the input (cf. H26ID and that have been used for computations with numerical data (cf., e.g., ||171 [TTD. 

In H17II a numerically stable version of the Buchberger-Moller algorithm has been derived which 
heavily relies on the usage of border bases in order to ensure numerical stability. An alternative algo- 
rithm based on border bases with similar stability behavior was investigated in [TJ. Both algorithms 
calculate bases of ideals of points which can be used to derive polynomial models from empirical 
datasets. 

It is well known that every Grobner basis with respect to a degree-compatible term ordering can be 
extended to a border basis (see 11201 p. 28 Iff]) but not every border basis is an extension of a Grobner 
basis. Moreover, not every order ideal supports a border basis even if it has the right cardinality. An 
example illustrating these two cases is presented in 11201 Example 6] . 

While the border basis algorithm in 1 12 Oil , which is a specification of Mourrain's generic algorithm 
[28], allows for computing border bases of zero-dimensional ideals for order ideals supported by a 
degree-compatible term ordering, it falls short to provide a border basis for more general order ideals: 
The computed border basis is supported by a reduced Grobner basis. The alternative algorithm pre- 
sented in [20, Proposition 5] which can potentially compute arbitrary border bases requires the a priori 
knowledge of the order ideal that might support a border basis so that the order ideal has to be guessed 
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in advance. As we cannot expect this prior knowledge, the algorithm does not solve the problem of 
characterizing those order ideals for which a border basis does exist. Further, as pointed out in 11201 p. 
284] , the basis transformation approach of this algorithm is unsatisfactory as it significantly relies on 
Grobner basis computations. 

It can be favorable though to be able to compute border bases for order ideals that do not necessarily 
stem from a degree-compatible term ordering. It has been an open question to characterize those order 
ideals of a given zero-dimensional ideal I that support a border basis. We will answer this question for 



degree-compatible order ideals by resorting to polyhedral combinatorics, mentioned in Subsection 1.2 
The restriction to degree-compatible order ideals is due to the design of the algorithm to proceed 
degree-wise and the authors believe that this restriction can be overcome as well. 

The preference of border bases over Grobner bases partly arises from the iterative generation of 
linear syzygies, inherent in the border basis algorithm, which allows for successively approximating the 
basis degree-by-degree. Moreover, the border basis algorithm is a linear algebra algorithm with a tiny 



exception, the final step which does not contribute to the inherent complexity though as its running 



time is polynomial in the input size. 

1.2. Polynomial method. In discrete mathematics and combinatorial optimization, polynomial sys- 
tems have been used to formulate combinatorial problems such as the graph coloring problem, the 
stable set problem, and the matching problem (see e.g., II13II for an extensive list of references). This 
well-known method, which Alon referred to as the polynomial method (cf. [[3l |4fl) recently regained 
strong interest. This emphasizes the alternative view of the border basis algorithm as a proof sys- 
tem which successively uncovers hidden information by making it explicit. In [14, Section 2.3] and 
11131 11211 infeasibility of certain combinatorial problems, e.g., 3-colorability of graphs is established 
using Hilbert's (complex) Nullstellensatz and the authors provide an algorithm NulLA to establish in- 



feasibility by using a linear relaxation. The core of the algorithm is identical to the L -stable span 



procedure used in the border basis algorithm, which intimately links both procedures. The difference 
is of a technical nature: whereas NulLA establishes infeasibility, the classical border basis algorithm as 
presented in [ 20 ] computes the actual border bases of the ideal. Another recent link between border 
bases and the Sherali-Adams closure (introduced in [33]), a method for convexification, was inves- 
tigated in II31II . The authors show that the Sherali-Adams procedure can be understood as a weaker 
version of the L -stable span procedure, which is an essential part of the border basis algorithm. As 
a consequence a new tighter relaxation than the Sherali-Adams closure is derived by exploiting the 
stronger L-stable span procedure. 

1.3. Applications of border bases. Surprisingly, it turns out that there are deep connections to other 
mathematical disciplines, and border bases represent the combinatorial structure of the ideal under 
consideration in a canonical way. Although the use of border basis as a concise framework is quite 
recent, introduced in ||191 [201 12111 and 11251 Section 6.4] , the concept of border basis is rather old and 
appeared in different fields of mathematics including computer algebra, discrete optimization, logic, 
and cryptography (independently) under different names. Originally, border bases were introduced in 
computer algebra as a generalization of Grobner bases to address numerical instabilities. Border bases 
have been successfully used since for solving zero-dimensional systems of polynomial equations (see, 
e-g-, U5l 1271 12811 ). which in particular include those with solutions in 0/1 and thus a large variety of 
combinatorial problems. 

Border bases have been also used to solve sparse quadratic systems of equations thus giving rise 
to applications in cryptography in a natural way. Such systems arise from crypto systems (such as 
AES, BES, HFE, DES, CTC variants, etc.) when rewriting the S-boxes as polynomial equations. The 
celebrated XL, XSL, MutantXL attacks, which are based on relinearization methods, are essentially 
equivalent to the reformulation-linearization-technique (RLT) of Sherali and Adams 1 13 3 1 1 and use a 
version of the Nullstellensatz to break ciphers. In fact the XL algorithm (see e.g., [22 , 8]) in its classical 
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form is actually identical to a level d Sherali-Adams closure of the associated system and therefore it 
is a border bases computation at its core. Motivated by the success of the aforementioned methods, 
border bases have also been used in crypto analysis and coding theory, see []6j • 

Another application of border bases is the modeling of dynamic systems from measured data (see 
e.g., IH71 1231 [TIP where the better numerical stability can be advantageous. 



Our contribution. A perennial problem in various applications is that the classical computation of a 
border basis depends on a degree-compatible term ordering and hence finds only the border bases sup- 
ported by an order ideal induced by the ordering. These border bases do represent only a rather small 
fraction of all possible border bases (see Example 2.6). Further, in this case the border basis contains a 
reduced Grobner basis, and for example the theoretically nice numerical/ symmetry properties are often 
lost, as the vector space basis of K [X] // (where 7 is a zero-dimensional ideal) might not be optimally 
suited for numerical computation. Moreover, for example when solving large systems of polynomial 
equations, it is desirable to guide the solution process by having a vector space basis that does actually 
have a nice interpretation attached. These bases often cannot be obtained by a degree-compatible term 
ordering. 

Techniques from commutative algebra, and in particular border bases (and alike) have been of great 
value to discrete mathematics and combinatorial optimization. This time we will proceed the other way 
around. We will apply the combinatorial optimization toolbox to the combinatorial structure of zero- 
dimensional ideals and we solve the aforementioned problems by freeing the computation of a border 
basis from the term ordering. Our contribution is the following: 

Polyhedral characterization of all border bases: We provide a complete, polyhedral character- 
ization of all degree-compatible border bases of any zero-dimensional ideal I: we associate 
an order ideal polytope P to I that characterizes all its degree-compatible border bases. The 
integral 0/1 points in the order ideal polytope are in one-to-one correspondence with degree- 
compatible order ideals supporting a border basis of I. (An order ideal is degree-compatible 
if it is a degree-wise complement of I, see Definition 2.3 ) This explicitly establishes the link 
between the combinatorial structure of the basis of the factor space and the structure of the 
ideal. Whether an order ideal supports a border basis is solely determined by the combinatorial 
structure of the order ideal polytope. A related result for Grobner bases of a vanishing ideal of 
generic points was established in H29II . 

Computation of border bases not coming from a term ordering: We outline an algorithm based 
on the classical border basis algorithm as defined in [20] and show how to compute bor- 
der bases for arbitrary degree-compatible order ideals without relying on a specific degree- 
compatible term ordering. Recall that not every order ideal supports a border basis and guess- 
ing the order ideal in advance is permissive. It is also not advised to search for an admissible 
order ideal with brute force as the combinatorial structure of the ideal can be so complicated 
that it is not even possible to easily guess a (non-term ordering induced) feasible order ideal, 
not even to mention to find one that has a preferred structure (see also the next point) . 

Finding preferred order ideals: Our approach is also able to compute order ideals maximizing a 
prespecified preference function (and compute their border bases). We will show that computing 
a preference-optimal order ideal supporting a border basis is NP-hard in general, which is 
surprising, as choosing the order ideal is merely a basis transformation. The NP-hardness does 
not come from the hardness of computing the L -stabilized span as the problem remains NP- 
hard in cases where the L-stabilized span is small enough to be determined efficiently. 

Computational feasibility: We provide computational tests that demonstrate the feasibility of 
our method. 

Applications: The structure of ideals and counting of border bases: Having the order ideal poly- 
tope available for a zero-dimensional ideal I, it is possible to examine the structure of the ideal 
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based on its border bases. One straight-forward application is counting the number of degree- 
compatible border bases for a zero-dimensional ideal I. 

Order ideals and determining those with maximum score do appear in a very natural way in combi- 
natorial optimization as the so called maximum weight closure problem (cf. e.g., [30]) and they have a 
variety of applications, e.g., in open-pit mining where any feasible production plan is indeed an order 
ideal; clearly we can only mine the lower levels after having mined the upper ones. A good survey as 
well as an introduction to the problem can be found in II18II . 

Other, more involved applications might arise, e.g., in computational biology where the structure of 
boolean network is inferred from the Grobner fan. Due to their higher numerical stability, border bases 
might be a better choice here and the order ideal polytope allows for a detailed study of the underlying 
combinatorial structure of the networks. 



Outline. We start with the necessary preliminaries in Section 2 and recall the classical border basis 



algorithm in Subsection 2.1 In Section 3 we introduce the order ideal polytope and establish the one- 



to-one correspondence between the 0/1 points of this polytope and degree-compatible border bases. 



We also derive an equivalent characterization that is better suited for actual computations. In Section 4 



we then show how the results from Section 3 can be used to compute admissible degree-compatible 
border bases using a preference function for the structure of the order ideal. We also establish the 
NP-hardness of the optimization variant of this problem. As an immediate application, we also outline 
how to count all admissible border bases of a zero-dimensional ideal. We conclude with computational 
results in lSection 61 and a few final remarks in lSection 7| 



2. Preliminaries 

We consider a polynomial ring K[X] over the field K with indeterminates X = {x 1; ...,x n }. For 
convenience we define x m := Yl}e[n] x j ' for m e N " and let T " := {* m I rn e N"} be the monoid 
of terms. For any d e N we let T" d := {x m e T n | \\m\\i < d} be the set of monomials of total 

degree at most d. For a polynomial p e K[K] with p = Xif=i a i xm ' we define the support of p to be 
supp(p) := {x mi | i e [I]} and similarly, for a set of polynomials P C K[X] we define the support of P 
to be supp(P) := {J peP supp(p). Given a term ordering o, the leading term LT^p) of a polynomial p is 
U a (p) := t with t e supp(p) such that for all t' e supp(P) \ {t} we have t > £J t'; the leading coefficient 
LC a (p) of p is the coefficient of LT CT (p). We drop the index a if the ordering is clear from the context. 
Recall that the degree of a polynomial p e K[X] is deg(p) := max x m esupp Q,) ||m||i. The leading form LF(p) 

of a polynomial p = $][=i a i xmi e K[X] is defined to be LF(p) = ||m-||=<J a i xTni where d = deg(p), 
i.e., we single out the part with maximum degree. Both LF and LT generalize to sets in the obvious way, 
i.e., for a set of polynomials P we define LF(P) := {LF(p) | p e P} and LT(P) := {LT(p) | p e P}. 

In the following we will frequently switch between considering polynomials M = {p\, ■ ■ ■ ,p s }, the 
generated ideal, and the generated vector space whose coordinates are indexed by the monomials in 
the support of M. We denote the ideal generated by M as (M) K [ X ] and the vector space generated by 
M as (M) K . For n e N we define [n] := {l,...,n}. All other notation is standard as to be found in 
H9l 12411 ; we have chosen the border basis specific notation to be similar to the one in II20II where the 
border basis algorithm was first introduced in its current form. 

Central to our discussion will be the notion of an order ideal: 

Definition 2.1. Let 6 be a finite subset of T". If for all t e C and t' e T n such that t' \ t we have 
t 1 e , i.e., is closed under factors, then we call an order ideal. Furthermore, the border 80 of a 
non-empty order ideal is the set of terms 80:= {xjt | ; e [n], t € 0} \ 0. As an exception, we set 
30 := {1} for the empty order ideal. 
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Recall that an ideal 7 c K[X] is zero-dimensional, if and only if K[X]/I is finite dimensional. The 
^■-border basis of a zero-dimensional ideal 7 is a special set of polynomials: 



Definition 2.2. Let = {t 1 ,...,t fX } be an order ideal with border 80 = {b 1 ,...,b v }. Further let 
^ = {g 1; ...,g v } C if[X] be a (finite) set of polynomials and let 7 C if[X] be a zero-dimensional ideal. 
Then the set ^ is an 0-border basis of 7 if: 

(1) the polynomials in 9? have the form g,- = b ; - — a if^i f° r 7 e M an< ^ a ij e ^> 

(2) m Km = i; 

(3) 7C[X] = 7 © (^) Jf as vector spaces. 

If there exists an ^-border basis of 7 then the order ideal supports a border basis of I. 

Note that the condition (^)jc[x] = 7 is a consequence of c 7, the particular form of the elements in 
c §, and K[X] = 7 © (0)k- More precisely, {^)k[x] + (^)jc is closed under multiplication by the x i} and 
hence it is an ideal. As it contains 1, it must be the whole ring, and hence {^)k[x] = J by the modular 
law. See | |21l Proposition 4.4.2] for another proof. 

In particular, an order ideal supports an ^-border basis of 7 if and only if K[X] = 7 © (0) K . 
Moreover, for any given order ideal and ideal 7 the ^-border basis of 7 is unique as bj has a unique 
representation in K[X] = I © {0) K for all j e [v]. Furthermore, as K[X] = I © {0) K it follows that 
\0\ = ciim{0) K is invariant for all choices of 0. The requirement for 7 being zero-dimensional is a 
consequence of this condition as well, as 80 (and in consequence 0) as part of the output of our 
computation should be finite. Clearly, as a vector space, 7 has a degree filtration, i.e., 7 = Uj eN 
where 7- 1 := {p e 7 | deg(p) < i}. For a set of monomials we define =l := {m e \ deg(m) = i} 
(and similarly for < instead of =). In the following we consider special types of order ideals, i.e., those 
that essentially preserve the filtration: 

Definition 2.3. Let 7 C K[X] be a zero-dimensional ideal and let c T" be an order ideal. We say 
that is degree-compatible (to 7) if 

dim (iT*) K = dim (T^- dim ^ 

for all i S N. 

Thus, the ^-border basis of a zero-dimensional ideal 7 with respect to any degree-compatible order 
ideal has a pre-determined size for each degree i e N. Intuitively, the degree-compatible order ideals 
are precisely those that correspond to degree-compatible orderings on the monomials. The important 
difference is that the orderings do not have to be term orderings. The definition above only requires 
local compatibility with multiplication as is an order ideal and thus downwardly closed and degree- 
compatible, i.e., if p,q are polynomials and deg(p) < deg(q) then p <q. The following example shows 



that the requirements of Definition 2.3 are not automatically satisfied by all order ideals 0. We also 
give an example of a degree-compatible ideal not coming from a term ordering. 

Recall that the order ideal associated to a Grobner basis of an ideal consists of all monomials not 
divisible by any leading term in the Grobner basis. 

Example 2.4 (A non-degree-compatible order ideal). With the degree lexicographic ordering on two 
variables, the set of polynomials 

3 2 2 3 

x l + x l x 2> x \ x 2> x l x 2> X 2 

is a Grobner basis of an ideal with associated order ideal consisting of all monomials of degree at most 
2: 

{1» x l, x 2> x \ > x l x 2> x 2,}- 
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The element + x±x 2 of the ideal enables us to replace x x x 2 with and thus obtain another order 
ideal of the ideal, which is not degree-compatible: 

{1, Xi, X 2 , X^, Xj, x^}. 

Example 2.5 (A degree-compatible order ideal not coming from a term ordering). The homogeneous 
ideal having a Grobner basis in the degree lexicographic term ordering 

has the associated order ideal 

{ 1 » x 1 , x 2 , x 1 x 2 , x 2 , x 2 } ■ 
Via a basis change, we obtain another order ideal 

{1, Xi, x 2 , x 1 ,x 2 ,x 2 }, 

which cannot come from a term ordering. The reason is that every Grobner basis of the ideal must 
contain a polynomial with degree-two leading term, which therefore must be a multiple of the first 
generator x^ + x x x 2 + x^ modulo higher-degree terms. The degree-two leading term cannot be XyX 2 
because if e.g., x 1 < x 2 then the term x\ is larger then x^x 2 . So the leading term must be either x^ or 
x^, which excludes the order ideal. 

Example 2.6 (Generic ideal). Let k and n be positive integers and let { a o he[n] je[/c] ^ e a lg e t ) raically 
independent real numbers over Q. Let I be the ideal of polynomials in the variables x-y, . . . ,x n which 
are zero on the points (a 1; , . . ., a n; ) for j e [fc]. Thus, the ideal is zero-dimensional, and i<C[X]/7 has 
dimension k. 

Every k distinct monomials form a complementary basis of I, since they are linearly independent on 
the k points (a 1; -, . . . , a nj ). An equivalent formulation of linear independence is that the determinant 
of the matrix formed by the values of the monomials on these points is non-zero. The determinant is 
indeed non-zero, as it is a non-trivial polynomial of the algebraic independent a i; - with integer coeffi- 
cients. 

In particular, every order ideal of size k is an order ideal of I. The degree-compatible order ideals 
are the ones where the monomials have the least possible degree, i.e., consisting of all monomials of 
degree less than I and k — ("t" 1 ) monomials of degree I, where I is the smallest non-negative integer 

satisfying k < ( n t ), i.e., there are at least k monomials of degree at most I. 

2.1. Computing border bases for term ordering induced 0. Without proofs, we recall the classical 
border basis algorithm introduced in [20] as it will serve as a basis for our algorithm. The interested 
reader is referred to | |21l 11911 for a general introduction to border bases and to II20II in particular for an 
introduction of the border basis algorithm. 

The classical border basis algorithm calculates border bases of zero-dimensional ideals with respect 
to an order ideal which is induced by a degree-compatible term ordering a. 

First, the border basis algorithm heavily relies on the following neighborhood extension: 

Definition 2.7. (cf. 11201 Definition 7.1 and paragraph preceding Proposition 13]) Let V be a vector 
space. We define the neighborhood extension of V to be 

y+ : =v + Vx 1 + --- + Vx n . 

For a finite set W of polynomials, its neighborhood extension is 

W + = W U Wxy U • • • U Wx n . 
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Note that for a given set of polynomials W such that (W) K = V we have {W + ) K = (W)i = V + as 
multiplication with x ; is a iC -vector space homomorphism. It thus suffices to perform the neighborhood 
extension on a set of generators W of V. 

Let F be a finite set of polynomials and let L c T" be an order ideal, then F n {L) K = {/ e F \ 
supp(/) c L}, i.e., F n (L) K contains only those polynomials that lie in the vector space generated by L. 
Clearly, for L = T" d we have {F) K n {L) K = (F)^ d . Using the neighborhood extension we can define: 

Definition 2.8. (cf. 11201 Definition 10]) Let L be an order ideal and let F be a finite set of polynomials 
such that supp(F) c L. The set F is L-stabilized if (F + ) K n (L) K = (F) K . The L-stable span F L of F is 
the smallest vector space G containing F satisfying G + n (L) K = G. 

A straightforward construction of the L-stable span of F is to inductively define the following in- 
creasing sequence of vector spaces: 

F :=(F) K and F k+1 := F+ n (L)* for k > 0. 

The union Ufc> F fc 1S tne L-stable span F L of F. 

The set L represents our computational universe and we will be in particular concerned with finite 
sets L c T". Note that L-stability is a property of the vector space and does not depend on the basis: 

Remark 2.9. The L-stable span of a finite set F depends only on the generated vector space (F) K , as 
(F + )k = (F) + k . 

In the following we will explain how the L-stable span can be computed explicitly for L = T" d . We 
will use a modified version of Gaussian elimination as a tool, which allows us to extend a given basis V 
with a set W as described in the following: 

Lemma 2.10. Q20l Lemma 12] Let V = {v 1; . . ., v r } c K[X] \ {0} be a finite set of polynomials such that 
LT(Vj) 7^ LT(Vj) whenever i, j e [r] with i ^ j and LC(v;) = 1 for all i e [r]. Further let G = {g 1; . . . , g s } 
be a finite set of polynomials. Then Algorithm 2.11 computes a finite set of polynomials W C K[X] with 
LC(w) = 1 for all w e W, LT(ux) ^ LT(u 2 ) /or any distinct u 1; u 2 S V U W, and (V U W) K = (V U G) K . 
(V, W may be empty.) 

Algorithm 2.11 (Gaussian Elimination for polynomials — GaussEl). 



Input: V, G as in Lemma 2.10 



Output: W c R[X] as in Lemma \2A0 

(1) Let H :=G and r? := 0. 

(2) IfH = then return W := {v r+1 , . . . , v r+r) } and stop. 

(3) Choose f S H and remove it from H. Let i := 1. 

(4) Iff = Qori>r + r\ then go to step 

(5) //LT(/) = LT(v;) then replace f with f - LC(/) • v £ . Set i := 1 and go to step @. 

(6) Set + Go to step @. 

(7) #7 7^ then put t] := rj + 1 and Zet v r+r) :=// LC(/). Go to step {^). 



We can now compute the L-stable span using the Gaussian elimination algorithm 2.11 



Lemma 2.12. 11201 Proposition 13] Let L = T" d and F C K[X] be a finite set of polynomials supported 
on L. Then Algorithm 2.13 computes a vector space basis V of F L with pairwise different leading terms. 

Algorithm 2.13 (L-stable span computation — LStabSpan). 
Input: F, L as in Lemma \2.12 



Output: V as in Lemma \2.12\ 

(1) V :=GaussEl(0,F). 

(2) W' := GaussEl(V, V+ \ V). 



(3) W:={weW' supp(w) c L} = {w e W' | deg(w) < d}. 

(4) Jf |W| > set V := V U W and go to step @. 

(5) Return V. 

The last ingredient that we need in order to formulate the border basis algorithm is the final reduc- 
tion algorithm. This algorithm basically transforms the output of the border basis algorithm to bring 
it into the desired form of a border basis by applying linear algebra steps only. It interreduces the 
elements so that they only have support in the leading term and 0. 

Lemma 2.14. 11201 Proposition 17] Let F be a system of generators of a zero-dimensional ideal I, let 
L be an order ideal, and let V be a vector space basis of F L with pairwise different leading terms and 



■— L \ LT( V) such that d c L. Then Algorithm 2.15 computes the -border basis <§ of I. 



It is easy to see that (L) K = F L © {G) K , so the algorithm will do a form of Gaussian elimination on 
the basis V to obtain a new basis V R with all non-leading terms supported on 0. Of course, this new 
basis will contain an ^-border basis. 

Algorithm 2.15 (Final Reduction Algorithm — FinalRed). 
Input: V, as in Lemma %2~14\ 
Output: <3 as in Lemma \2.14\ 



(1) Let V R := 0. 

(2) IfV = then go to step @. 

(3) Let veV such that v has minimal leading term. Put V := V \ {v}. 

(4) Let H := supp(v) \ (LT(v) U 0). 

(5) If H = then append v/LC(v) to V R and go to step @. 

(6) For each h e H choose w h e V R and c h ^K such that LT(w h ) = h and h ^ supp(v — c h w h ). 

(7) Set v '. — v — X^eH c h w iv append v/ LC(v) to V^, and go to step (j^J). 

(8) Return <g := {v e V R \ LT(v) e 30}. 

We will now formulate the border basis algorithm. 
Proposition 2.16. 11201 Proposition 18] Let F c K[X] be a finite set of polynomials that generates a 



zero-dimensional ideal I = (L^rxi- Then Algorithm 2.17 computes the 0-border basis of I. 



Algorithm 2.17 (Border basis algorithm — BBasis). 



Input: F as in Proposition 2.16 
Output: ^ as in Proposition 2.16 

(1) Let d := max^ eF deg(/). 

(2) V = {v 1; ...,v r } :=LStabSpan(F,T^). 

(3) Let := T^ d \ {LT( Vl ), • • • , LT(v r )}. " 

(4) Ifd0£ then set d := d + 1 and go to step ||). 

(5) Return <g := FinalRed(V, 0). 



It is worthwhile to note that step Q in Algorithm 2.17 is essentially for testing if the L-stable span 
is large enough to support an ^-border basis. 

The rationale for searching such a large span is due to the following proposition that serves as a 



stopping criterion. It is obvious from its Corollary 2.19 that the span is indeed minimal in the sense 



that it is the smallest span that contains a degree-compatible order ideal that supports a border basis 
(and thus all such degree-compatible order ideals) . 

Proposition 2.18. 11201 Proposition 16] Let L be an order ideal. Further let I be an L-stabilized generating 
vector subspace of a zero-dimensional ideal I Q K[X], i.e., I + n (L) K = I and ^ R = I. If is an order 

ideal such that {L) K = / © {0) K and 80 Q L then supports a border basis of I. 



We obtain the following corollary: 



Corollary 2.19. Let I be an T n <d -stabilized vector space satisfying I + (t"^) = (T" d ) . Then 



Proof. We apply Proposition 2.18 with the choice L := T^, I := (i) K[X] and := \LT(7) where the 
leading terms are with respect to any degree-compatible term ordering. Clearly, {^<^ K = I © (@)k- 
The condition 7+ = {^<d) K ensures that 6 consists of monomials of degree less than 

d, so 80 c T" Hence the proposition applies, and we obtain 7C[X] = 7 © (0) K . Together with 
{^k d ) K = 1® Wk this gives 7n (T^ d ) K = 7~. □ 

The border basis algorithm decomposes into two components. The first one is the calculation of the 
L -stable span for L = T" d for sufficiently large d (this is the main 'work') and the second component 



is the extraction of a border basis via the final reduction algorithm Effectively, after using any degree 



compatible term ordering in order to compute the L -stable span one can choose a different ordering 
with respect to which the basis can be transformed. This approach is very much in spirit of the FGLM 
algorithm, Steinitz's exchange lemma, and related basis transformation procedures. In Remark 2.9 it 
is shown that L -stability does not depend on the basis and thus any sensible basis transformation that 
results in an admissible order ideal is allowed. In the following we will characterize all admissible 
order ideals. 



3. The order ideal polytope 

We will now introduce the order ideal polytope P (a 0/1 polytope) that characterizes all order ideals 
that support a border basis (for a given zero-dimensional ideal 7) in an abstract fashion independent 



of particular vector space bases for the stable span approximation. Its role will be crucial for the later 



computation (Algorithm 4.1) of border bases for general degree-compatible order ideals. In the first 
subsection, we introduce the polytope in an abstract, invariant way highlighting its main property 
that its integral points are in bijection with degree-compatible order ideals supporting a border basis 



(Theorem 3.2). In the second subsection, we give a more direct reformulation targeted to actual 



computations. 



Given d G N in advance for which Algorithm 2.17 stops, LStabSpan computes the actual border 



basis (and some unused extra polynomials). The computation of the border basis is performed by 
Gaussian elimination applied to the matrices obtained within LStabSpan. The order ideal for which 
we effectively compute the border basis is solely determined by the pivoting rule when doing the 
elimination step. Degree-compatible term orderings ensure that we obtain an order ideal for which a 
border basis exists (given the correct d). If we would now pivot arbitrarily, which effectively means 
permuting columns, it is not clear that, first, the resulting ideal is an order ideal and, second, that it does 
actually support a border basis. If we remain in the setting of degree-compatible order ideals (which 
means that we only permute columns of monomials with same degree) then we face the combinatorial 
problem that when permuting a column such that it will end up being an element of we also have 
to ensure that all its divisors will also end up in 0. On the other hand, due to the degree-compatibility 



constraint (see Definition 2.3 1 we must choose an exactly determined amount of elements for each 
degree that will end up in 0. We will show now how to transform this combinatorial problem of 
faithful pivoting into a polyhedral setting. We obtain a 0/1 polytope P, the order ideal polytope that 
characterizes all admissible degree-compatible order ideals that support a border basis for the ideal at 
hand. 



3.1. Theoretical point of view. 
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(3.1) 
(3.2) 

(3.3) 
(3.4) 



z mi >z m2 Vm lj m 2 <ET n : m 1 \m 2 
^ z m = dim(Tl i ) K -dimI^l^ i - 1 Vi 



meT" . 



J] z m < dim ([/ U I^/l^- 1 ) R - dim 1^ jl^' 1 

Vi, [7 c : |[7| = dim (r n =i ) K - dim I^/l- 1 ' 1 
z m e[0,l] Vmef 



met/ 



Figure 3.1. Order ideal polytope P(7) 



Definition 3.1. Let 7. be a zero-dimensional ideal. Its order ideal polytope P(7) is given by the system of 



inequalities in |Figure 3.1| 

The order ideal polytope is actually a finite dimensional polytope as all the z m ar e when the d egree 



gives 



of m is large enough. Indeed, for large i, we have (T-')^ = 7-'/7-' _1 and hence Condition 3.7 
z m = for every m of degree i. 

We are ready to relate the order ideal polytope with order ideals. From now on, let A(7) denote the 
set of degree-compatible order ideals of a zero-dimensional ideal 7. 

Theorem 3.2. Let I be a zero-dimensional ideal. There is a bijection between the set A(7) of its degree- 
compatible order ideals and the set of integral points of the order ideal polytope of I. The bijection is given 
by 

£: zeP(/)nf" ^ 0{z) := {meT" | z m = 1}. 
Proof In fact, we will see that the order ideal polytope is defined exactly to this end. 



First, the integral solutions z of Condition 3.4 are exactly the 0/1 points, i.e., the characteristic 
vectors of sets of terms 0{z) := {m G T" | z m = 1}. 



Second, it is easy to see that Condition 3.6 means that 0{z) is indeed an order ideal, as whenever 



mi | m 2 and m 2 e 0{z), i.e., z m2 = 1, then it follows that z mj = 1, i.e., m 1 e 0{z) as well. 

In the third step we will provide an algebraic characterization of Conditions 3.2 and 3.3 Clearly, 
ICondition 3.21 can be rewritten to 



\0{z) =i \ = dim(T^) K -dim7- i /7- ; - 1 



We will now show that Condition 3.3 is equivalent to 
(3.5) {0{zT i ) K (M- i jl- i - 1 ={0}, 

i.e., the image of 0{z) =l is linearly independent in the factor ( T =;) K /(- f "7^~ !_1 ) ■ 
Similarly as above, [Condition 3.3| can be rewritten to 

| U n 0{z) =i | < dim (U U I- 1 1 1- 1 ' 1 ) R - dim I^jl- 1 ' 1 , 

i.e., the size of U n 0(z)~ l is at most the dimension of the vector space generated by the image of U 
in the factor (T^^^j ^7- ! j 7- !_1 ) . This is obviously necessary for the image of 0(z)~ l to be linearly 
independent in the factor. (Here the size of U does not matter.) 

For sufficiency choose U := 0(z)~ l . Then the dimension of the vector space generated by the image 
of 0(z)~ l is at least \0(z)~ l \, so the image of 0(z)~ l is independent. 
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So far we have proved that the integral points of the order ideal polytope correspond bijectively to 
order ideals with the properties 

\0{z) =i \ = dim (T n =i ) K - dim I^jl-^ 1 

and 

(^r^n^/i^- 1 = {0} 

for all i. 

Lastly, we will show now that this is equivalent to |^(z) =! | = dim^T"^ — dim and 
I © {6{z)) K = (T n ) K and thus the assertion follows. 

By dimensionality it follows, \0(z) =i \ = dim (t^) - dim I^/l- 1 ' 1 and (^f^n^'/j-'" 1 = 
{0} for all i together are equivalent to © {^0{z) =l ^ = (T™J for all i. For brevity we will 

omit the phrase 'for all V. Using the filtration argument, the latter is equivalent to I- 1 © ^&(z)- l ^ R = 
('" , <i)jf" ^ a a dimension argument on embeddings, this is further equivalent to 7 © {0{z~)) K = (T n ) K 
and dim/- ! + dim (&(z)- l \ R = dim ^T< { \ . Finally, filtrating the dimension by degree shows that the 
latter is equivalent to J © (0(z)) K = (T n ) K and dim J-'//"' -1 + dim {^0{z) =i ) K = dim (T^) . □ 

3.2. Computational point of view. Throughout this subsection we assume that M = IJieN — ^ M 
is a finite set of polynomials with degree-filtration {M i \ i e N} that generates a zero-dimensional 
ideal (M) K r x i C K[K] such that all p e M ; have degree i. Furthermore, for each i e N we have an 
enumeration M; = {p i; - | j e [fc;]} with k t e N. As seen in Subsection 2.1 an important component 



of the border basis algorithm is the computation of L-stable spans with respect to some computational 
universe L c T n . Computing a border basis with respect to a different order ideal is merely a 
basis transformation of the vector space obtained from the LStabSpan procedure. In the following 



we assume that L is of the form L = T" d where d e N is such that Algorithm 2.17 stops. In view of 
Remark 2.9 and Definition 2.3 we have that for a given computational universe L either all border bases 
(supported by degree-compatible order ideals) are contained in L or none. Furthermore we assume that 
M is L -stabilized and has in particular a convenient form. Effectively one might want to think of M 
being the output of the LStabSpan procedure, which is then brought into the following reduced form: 

Definition 3.3. Let M be a finite set of polynomials of degree at most I for some I e N. Then M is 
in canonical form if the leading term of any element of M does not occur in the other elements. Here 
we can freely choose leading terms of the polynomials with the only constraint that they have to be 
maximal-degree terms. 

We give a visual interpretation of the definition. The coefficient matrix A e K <t of M is the 
matrix where the rows are the elements of M, the columns are all the monomials of degree at most I, 
and the entries are the coefficients of the terms in the elements of M. We use the convention that for 
terms t x , t 2 with deg(tj) > deg(t 2 ) we put column t x to the left of column t 2 . Similarly, we put leading 
terms of a polynomial to the left of the other terms. Now M is in canonical form if the matrix A has the 



structure as depicted in Figure 3.2 i.e., it consists of degree blocks and each degree block is maximally 
inteneduced. 

The degree blocks correspond to the leading forms of the polynomials in M. Any finite set can be 
brought into canonical form by applying Gaussian elimination and column permutations of terms with 
same degree if necessary. In particular, the output of the LStabSpan procedure can be easily brought 
into this form. The following lemma summarizes the basic properties of a set M in canonical form: 

Lemma 3.4. Let M be in canonical form and L-stabilized with L = T" rf . Let T" d C (M) K + T" d _ r Then 
the following holdfor all i e [d]: 

(1) (M)|| x] /(M)| ! - ] 1 = (LF(M i )) Jf 

li 



A = 



1 
1 


-k 























1 
1 













































1 
1 


* J 



Figure 3.2. canonical form 



(2) (M)|f x] = (u<iM ; ) K 

(3) (M;)k = and thus (M t )f C ((J 



^0<;<i-l M ; 

Proof. We first show that (M;)^' = for all i e [d]. Let i e [d] be arbitrary and observe that each 
nonzero element p e M ; has degree i. As M is in canonical form, the polynomials in M; are interreduced 



(see the matrix in Figure [3^2] for Definition 3.3) and thus we also obtain each nonzero element p e 
{M i ) K has degree i. 
By Corollary 



2.19 



(M) K[X] n (L) K = (M) K . Hence (M)| 



(M)| ! for i S [d]. Now the statements 

□ 



of the lemma are obvious consequences of M being in canonical form. 

The following lemma provides us a practical way to compute the sizes of the degree components of 
degree-compatible order ideals, which are the same for all order ideals of a given ideal. 

Lemma 3.5. Let M be in canonical form and L-stabilized with L = T" d and d = max me3e deg(m). Let 
T" d c (M} K +T< d _ r Further let be an order ideal of (M) K [x] . Then is degree-compatible if and only 

if _ _ 

\0~ l \= dim (l~ 1 ) k - dim (LF(Mj)) K 

/or every i G [d]. 

r <i-i ^ ^ LF ( M .)^ by Lemma 



Proo/ In view of Definition 
where 7 := (MWi. 



2.3 



it suffices to observe that 7- 1 j I 1 



3.4 



(1) 

□ 



We are ready to provide a reformulation of the definition of order ideal polytopes, which is better 
suited for actual computations, partly as it only involves direct matrix operations via replacing dimen- 
sions with ranks of subsets: 

Lemma 3.6. Let M be L-stabilized and in canonical form with L = T n <d and (M)| d j ' {M)^ d ~ l = (T" d ) . 
Then the order ideal polytope P(M, L) c [0, 1] L of (M) K pn is given by the system of inequalities in 



Figure 3.3 



Proof. Let I := (M) K ^. We successively transform the defining inequalities of the order ideal polytope 

Figure 3.3 The reformulation is mostly based on I- 1 jl- 1-1 = 



m 



Figure 3.1 



into the desired form of 



(LF(M ; )) K from Lemma |3T4][lj). 

r <d-l 



First, as I- d jl'- 



f^=d)jc' we can remove tne variables z m with m degree at least d together 
with the inequalities involving them. These variables are always zero. 
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(3.6) 
(3.7) 

(3.8) 



z m 1 


>z m2 Vm 1) m 2 eL:m 1 |m 2 






= dim(L ! ) Jf -dim(LF(M i )) Jf 


Vi e [d - 1] 


meL=< 








> |[7|-rk(t7) 




meL/ 


Vie [d-l],[/CL =, ' : |[/| = 


dim (LF(M ; )) k 


z m 


e[0,l] VraeL 





Figure 3.3. Order ideal polytope P(M, L). In Condition 3.8 the matrix U is the induced 



sub-matrix of LF(M ; ) with column monomials only in U. 



Second, we replace all occurrences of I- 1 jl- 1 1 with (LF(M;)) K (or simply LF(M;)). This almost 



results in the inequality system of |Figure 3.3| with the only difference that instead of |Condition 3.8 
have 



we 



(3.9) Yt z ™ <dim([/ / ULF(M i )) K -dim(LF(M i )) K 

Vi e [d - 1], U c L =i : \U'\ = dim (L =i ) K - dim (LF(M;)) K , 

where we have deliberately replaced U with U'. 

We will show that the difference of Conditions |3 . 7| and |3 .8| is equal to Condition 3.9| with the choice 
U' := L =i \ U, which has size \U'\ = dim (l =i ) - dim (LFfM;))^ . This will finish the proof. 

Let [/ c L =I as above and compute the difference of Condition 3.7 and Condition 3.8 We obtain 



z m< dim (L~ i ) K - dim (LF(M ; )) K - \U\ + rk(t>). 

meL =i \£/ 



It is easy to see that \U\ = dim 



rk(L/) = dim 



. We claim that it suffices to show that 
LF(M i )u(L= i \U)) K 



Indeed, using this we can rewrite the inequality as 



^ / -x ( L=i )r (hF(M i )Ll( < L- i \U)) l 

2^ z m < dim {L~ l ) K - dim (LF(M;)) K - dim , , _,. , .* + dim ■ 



meL= l \U 



dim (LF(Mi) U (L =i \ U)) R - dim (LFfM;))^ , 



{L =i \U) K 



which is d3T9] ) for U' := L _1 \ L7 as claimed. 

We will show now that rk(t/) = dim (lF(M ; ) U (L = ' \ C/)) K - dim (l =i \ Let B denote the 

matrix obtained when writing the elements in LF(M ; ) as rows and let L =l index the columns. Clearly, 
(i =! ) K = (Bi<=L =i as a vector space and (U') K = (6[J / iCe f as a sub vector space with U' := L =, \U. 

We obtain (l =i ) K j{U') K = ®i eL ^\u' Ke t- Note that { L KMi)) K £ (L =i ) K and thus we can consider 

(LF<iM i )UU / ) K l(U / ) K C (L =i ) K l(U') K . Now dim ((LF(M i )U[/ / ) K /( C/Ok) = rk (C>) where [/ is 
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obtained from B by removing the columns in U' . We obtain 

(LF(M;) U U') k 



rk([7) = dim 
= dim 



(h¥(M i )U(L= i \U)) K 



(L =i \U) K 

and thus the result follows. □ 
4. Computing border bases using the order ideal polytope 



In the following we explain how Theorem 3.2 can be used to actually compute border bases for 



general degree-compatible order ideals. We cannot expect to be able to compute a border basis for any 
degree-compatible order ideal, simply as such a basis does not necessarily exist. Having the order ideal 
polytope at hand something slightly more subtle can be done: By choosing a linear objective function 
c € Z T and optimizing it over the order ideal polytope, we can actually search for an order ideal with 



preferred monomials in its support (see Subsection 4.2). Having the order ideal polytope available we 



can also count the number of degree-compatible border bases that exist for a specific ideal. Before we 
can address this application though, we will first show how to obtain an ^-border basis for e A(J) 
where I C K[X] is a zero-dimensional ideal. 

4.1. Computing border bases for e A(7). As the computation of the L -stable span of a set of 



generators M is independent of the actual chosen vector space basis (see Remark 2.91, we can adapt 



the classical border basis algorithm (Algorithm 2.17) to compute border bases for general degree- 
compatible order ideal. We first determine the right computational universe L = T" , for some d e N 
such that the associated L -stable span M contains all border bases. In a second step we optimize over 
the order ideal polytope P(M,L) and then perform the corresponding basis transformation. We will 
first formulate the generalized border basis algorithm by adding two steps after (|3]) in Algorithm 2.17 to 



the classical border basis algorithm, formulate the missing parts, and then prove its correctness: 

Algorithm 4.1 (Generalized border basis algorithm — BBasis). 

Input: F a finite generating set of a zero-dimensional ideal. 
Output: ^ a border basis of the ideal. 

(1) Let d := max /eF {deg(/)} and put L := T^ rf . 

(2) V = {v 1; . . . , v r } := LStabSpan(F, L). 

(3) IfT n =d £ LT(V) then setd:=d + l and put L := T^ d and go to step 

(4) Choose e A((V) K[X] ) (oze P(V, L) n Z L and = 0{z)). 

(5) Let ^ := BasisTransf ormation(V, 0). 

Note that step ([3]) is a convenient way to quickly check whether V is already L -stabilized and if 
it contains all degree-compatible order ideals. In this augmented algorithm we added the steps ([4]) 
and The first step will be extensively discussed in Subsection 4.2 as there are various ways to 



determine e A((V) K [ X ]) and this is precisely one of the main features, i.e., to choose the order ideal 
more freely. Note that by Theorem |3.2| we already know that does support an ^-border basis of 
(^)k[x] ( as (F)k [x] = (^)k[x]) an d our task is now to actually extract this basis from V. This extraction 
is performed in step ([5]). 

Let A be a matrix representing a set of polynomials M where the columns correspond to the mono- 
mials in some fixed ordering. Let the head (short: Head(a)) of a row a of A be the left-most monomial 
in the matrix representation whose coefficient is non-zero. Note that the notion of head replaces the 
notion of leading term of a polynomial as we do not (necessarily) have a term ordering anymore. The 
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main idea is to reorder the columns of V and then to bring V into a reduced row echelon form such 
that no m e is head of a row of the resulting matrix — a classical basis transformation: 



Lemma 4.2. Let L = T™ . with <eN, let V be a finite set of polynomials satisfying (V) K = (V}jcpn n (L) K 



returns 



and let = {t 1 ,...,t fX } be an order ideal with 30 Q L and e A((V) JC [ X ]). Then Algorithm 4.3 
the -border basis <g of {V) K [x\- 

Proof. First, the algorithm finds I from M. As e A((V) K -m) we have that supports a border basis 
and in particular we have K[X] = (V) K [ X ~] © (&)k> an d hence (L) K = (V) K © (0) K by the modular law. 
We will now show that Condition (1)| of Definition 2^2 is satisfied. This in turn follows from the fact 
that the algorithm creates every element gj of <g to have the form 



Sj = b j -Y i a iJ t i 



i=l 



with a t j e K. By construction b ; - e 30 and thus ^ is an ^-border basis of (V) 
Algorithm 4.3 (Basis transformation algorithm — BasisTransf ormation). 



□ 



Input: V, as in Lemma 4.2 



Output: <3 as in Lemma \4.2\ 

(1) Set I ■— max meM deg(m). 

(2) Permute the columns of V such that t is right ofm in the matrix representation of V for all m 
and t e 0. 



T" 



(3) Reduce V using Gaussian elmination: like Algorithm 2.11 but use Head instead of LT and the 
coefficient of the head instead of LC. Let <£' be the result. 

(4) Let <S := {g e <§' : Head(g) e 30}. 

(5) Return C S. 



Note that Algorithm 4.3 works for any order ideal that supports a border basis of {V) K ^, i.e., 
also those that are not necessarily degree-compatible. When the order ideal is known to be degree- 
compatible, it is enough to do the permutations in each degree block in the first step, and then use 



Algorithm 2.11 in the second step. 



We will show now that Algorithm 4.1 computes an ^-border basis for e A(7). 



Proposition 4.4. Let F = {/i,...,/ s } Q K[X] be a finite set of polynomials that generates a zero- 
dimensional ideal I = Then Algorithm 4.1 computes the -border basis of I for any (chosen) 



Proof. Whenever we reach step Q in Algorithm |4.1| we have that V is L -stabilized for some L = 
N and it contains all degree-compatible order ideals supporting a border basis, i.e., all 
A(7). Observe that J = (F) K [x] 

l ^ ^ T „^ 



T^ d with d 



an ^-border basis of (F) K 
L 



Note that step (3|) ensures 3 



{ 



Ik 



and thus, by Lemma 4.2 it follows that ^ is indeed 

C L via {V)f/{V)f- 



□ 



An improved version of the border basis algorithm has been also considered in [20]. Basically, the 
improvement can be traced back to considering more restricted computational universes L that arise 
from choosing L to be the smallest order ideal that contains the support of the initial system and 
then successively extending it using the + operation. This improvement due to the restriction of the 
computational universe L cannot work in our setting anymore: Suppose that V is L -stabilized with 
respect to some computational universe L T" d for all d e N (i.e., L is not obtained by bounding the 
total degree of the monomials in T) and contains the order ideal that is induced by the chosen degree- 
compatible term ordering in the classical border basis algorithm. Then the associated polytope P(V,L) 
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would only contain a subset of all possible degree-compatible order ideals, as we might be lacking 
monomials that we need to represent certain alternative choices of 0. If a subset of all admissible 
degree-compatible order ideals is sufficient, then the same optimizations can be applied though. 



4.2. Computing preferred border bases. Let V and L be as obtained after step ^ in Algorithm 4J. 



As shown in Theorem 3.2 and Lemma 3.6 the order ideal polytope P(V, L) characterizes all degree- 
compatible order ideals that support a border basis of {V) K [ X y Every z e P(V, L) n Z L induces an order 
ideal 0{z) which supports an ^(z)-border basis of (V) K [ X y This also shows that we cannot expect 
that every order ideal supports an ^-border basis of (V)k[x] as tne characterization is one-to-one. 
The natural question is therefore how to specify which order ideal should be computed, i.e., which 
z e P(V, L) n Z to choose. As we cannot always get what we would like to have, it suggests itself to 
specify a preference, i.e., which monomials we would like to be contained in 0{z) and which ones we 
would rather not. As the coordinates of z are in direct correspondence with the monomials in L we can 
define: 

Definition 4.5. A preference is a vector c e Z 1 which assigns a weight to each monomial m e L. If 
z G P(V, L) n Z L , then cz is the score or weight of z. 

As P(V, L) c [0, 1]" is a polytope we can optimize over P(V, L) n Z 1 and compute an element z e 
P(V, L) n Z L that has maximal score, i.e., we can compute z e P(V, L) n Z L such that 

cz Q = max{cz | z e P(V, L) n Z L }. 

In this sense a preference is an indirect way of specifying an order ideal. For certain choices of c € Z 
though it can be hard to compute an order ideal that maximizes the score as we will show now. 



5. Complexity of finding preferred order ideals 

In this section, we show that finding a weight optimal, border basis supporting order ideal of a zero- 
dimensional ideal given by generators is NP-hard (Theorem |5.3[ ). Note that this also translates to large 
classes of other choice functions as P(V,L)nZ L is a 0/1 polytope and thus its extremal points are given 
by an inequality description. So any choice function that implicitly asks for a fc-clique (which we will 
use in our reductions) can be replaced by the corresponding linear function and hardness also follows 
in this case. The hardness for computing a weight optimal order ideal is unexpected in the sense that 
we merely ask for a basis transformation. On the other hand it highlights the crucial role of order ideals 
in describing the combinatorial structure of the ideal. 

As an immediate consequence it follows that it is rather unlikely that we can obtain a good charac- 
terization of the integral hull conv(P(V, L) n Z 1 ) and we will not be able to compute degree-compatible 
order ideals that support a border basis and have maximum score efficiently unless P = NP. This shows 
that not only computing the necessary liftings of the initial set of polynomials via the LStabSpan pro- 
cedure is hard but also actually determining an optimal choice of an order ideal once an L -stable span 
has been computed. As mentioned before, this is in some sense surprising as the actual interference 
has been already performed at that point and we are only concerned with choosing a nice basis. From 
a practical point of view this is not too problematic as, although NP-hard, computing a weight optimal 
order ideal is no harder than actually computing the LStabSpan in general. For bounds on the degree 
deN needed to compute border bases see, e.g., 11131 Lemma 2.4]; the border basis algorithm gener- 
ates the Nullstellensatz certificates and is therefore subject to the same bounds. Further, state-of-the-art 
mixed integer programming solvers that can solve the optimization problem such as scip ([]2]]), cplex 
( II10II "). or gurobi ( 01611 ) can handle instance sizes far beyond the point for which the actual border 
bases can be computed. Very good solutions can also be generated using simple local search schemes 
starting from a feasible order ideal derived from a degree-compatible term ordering. 
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Figure 5.1. Order ideal computation as minimum cut problem. Capacities denoted 
next to the arc and weight c u denoted next to the respective node u. The dashed arcs 
have capacity oo. 



5.1. Fast without constraint. Determining an order ideal of maximum score in a computational 
universe L without having any additional constraints can be done in time polynomial in |L|, as was 
shown in II30II . One simply transforms it into a minimum cut problem in graphs (see, e.g., [32]): Let 
c e Z L be a preference vector. Define a graph F := (V, E) with V := L U {s, t} and E := {(iz, v) | u, v e 
L and v | u}, i.e., whenever v | u we add an arc from u to v. In fact, it is enough to have an arc when 
u = vx for some variable x, i.e., to consider the transitive reduction of E. Define 

E := E U {(s, u)|u£L,c„>0}U {(u, t) | u s= L, c u < 0}. 

Further we set all the capacities of the arcs with both vertices in L to oo, for any arc (s, u) with ueL 
we set the capacity to c u , and for any arc (iz, t) with iz e L we set the capacity to \c u \; let zc(iz, v) denote 



the capacity of the arc (iz, v) e E. An example is depicted in Figure 5.1 

For U,W c V, we define C(L/, W) := E( UjW )e[/xw w )- A fotj-azt (S,S) is a partition S U S = V 
of the vertices of V with seS and t£S and the weight of the cut is C(S,S). 

We would like to compute an order ideal contained in L with maximum score: 



max 



u&0 



c L order ideal 



Observe that, if (S,S) is a cut in T of finite weight, there exist no arc (iz, v) e £ with ueS and v e S. 
Therefore, (S,S) is a cut in T of finite weight if and only if S \ {s} is an order ideal. We can therefore 
rewrite the optimization problem as follows: 



max 



c L order ideal j = max{C({s}, 0) - C{0,{t}) \ c L order ideal} 

= max{C({s}, L) - C({s}, L\0) — C{0, {t}) \ c L order ideal} 
= C({s}, L) - min{C({s}, L\0) + C{0,{t}) \0QL order ideal} 
= C({s}, L) - min{C({5} U 0,(L\0)U {t}) \0QL}. 

The last line asks for a minimum weight cut in the graph F. Note that we can indeed drop the condition 
that has to be an order ideal as it is guaranteed implicitly by all finite weight cuts as explained above. 
It is well-known that the computation of a minimum cut can be performed in time polynomial in the 
number of vertices and arcs (see [34]). Thus we can indeed compute an order ideal of maximum 
score efficiently in this case. 
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5.2. NP-hard with constraints. So far we did not include the additional requirements as specified 



by the order ideal polytope (see Figure 3.3), in order to obtain degree-compatible order ideals that do 
actually support a border basis of the ideal I under consideration. Unfortunately, when including these 
additional requirements, the problem of computing an order ideal of maximum score is NP-hard as we 
will show in the following. 



We will show NP-hardness by a reduction from the MAX-CLiQUE| problem, which is well known to be 



NP-complete (see, e.g., II15II or Hill GT22]). Given a graph r = (V, £)> recall that a clique C is a subset 
of V such that for all distinct u,veC we have (u, v) € E. 

Max-Clique: Let r = (V, E) be a graph. Determine the maximum size of a clique C contained in F. 

We will in particular use the following variant: 

/c-Clique: Let r = (V, E) be a graph. Determine whether r contains a clique C of size k. 

Note that if r contains a clique of size k for some k e N, then so does it for any 1 < I < k. Further 
the maximum size of a clique is bounded by |V|. It follows that already testing whether F contains a 



clique C of size k for some given k e N has to be NP-complete, as otherwise, we could solve the Max 



Clique problem with 0(log 2 |V|) calls of the algorithm that solves the test variant via binary search. In 
[18, Discussion after Definition 3.2] it was indicated that determining a maximum weight order ideal 



of a pre-defined size is NP-hard by a reduction from |Max-Clique| While already indicating hardness 
due to one additional cardinality constraint, this is a slightly different problem than the one that we 
are facing here: in addition to the size constraint, the order ideal has to be degree-compatible. Thus 
we have a cardinality/degree-constraint for each degree of the monomials (Constraints ( |3.7| ) and ( |3.9| ) 
in Figure 3.3). Further, the degree constraints are not completely independent of each other. So it 
could be a priori perfectly possible that this particular restriction of the problem can be actually solved 
efficiently, which is not the case as we will see soon. 
We consider the following optimization problem: 

Max-Bounded order ideal: Let L = T" d for some d e N and c e Z L be a preference. Further let d t e N 
for i e [d — 1]. Determine an order ideal O c L with \@q 1 \ = d ; for all i e [d — 1] and maximum score, 
i.e., 

c(? _1 (^ )) = max{c(rV)) \0QL order ideal and \0 =i \ = d f for all i e [d - 1]}. 



By a reduction from /c-Clique we obtain: 
Theorem 5.1. IMax-Bqunded order idealI ls NP-hard. 

Proof. Let r = (V,£) be an arbitrary simple graph and choose k e [|V|]. We consider the polynomial 
ring K[x v | v S V]. Let I = T^. Further choose d 1 — k and d 2 = k(k + l)/2. Define 

{1, if m = x u x v and (iz, v) e E; 
0, otherwise. 

A vertex v e V is represented by the degree-one monomial x v and each edge (iz, v) e E is represented 
by the degree-two monomial x u x v . In some sense we extended the graph (V, E) to the complete graph 
K"| V | that has |V|(|V| — l)/2 edges as we consider all possible degree-two monomials that correspond to 
the edges; those edges that do not belong to E have weight though. 

The order ideals satisfying the size constraints are those consisting of the monomials 1, x v ., x v .x v . 
for 1 < i,j < k for some distinct vertices v lv . . ,v^. The score of the order ideal is the number of edges 
in the subgraph spanned by the x v . 

Thus when maximizing c we ask for a k -vertex subgraph with maximal number of edges. The 
subgraph is given by the degree-one monomials contained in the order ideal, which we denote by , 
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i.e., 

C := {v e V | x v e 

So the maximum score is equal to — l)/2 if and only if r contains a clique of size k. Thus we can 
solve /c-Clique efficiently if we can solve Max-Bounded order ideal efficiently and so the latter has to 
be NP-hard. □ 

Finally, it remains to show that for every graph r = (V, E) and k e [\V\] there exists a system of 
polynomials Fw \,k — K[x v \ v e V] spanning a zero-dimensional ideal such that solving the 



IDEAL OF IDEAL 



Max-Order 



problem for F\v\,k solves the /c-Clique problem for T. For this, we construct an ideal 
encoding all fc-cliques of the complete graph on n vertices: Let n e N and k e [n] and define 

ie[n] 

F n ,k :={v J -|;e[n-fc]}UT^ 3 . 

We consider the ideal generated by F nk . We show that its order ideals are in one-to-one correspondence 
with the fc-element subsets of the set of n variables x ls . . . ,x n as stated in the following lemma. 

Lemma 5.2. Let n e N and k e [n]. Then F nk generates a zero-dimensional ideal such that e 
A {( F n,k) K[x] ) if and only if =1 C with \0 =l \ = k, 0= 2 = {xy \ x,y S 6=\ and 6 =l = 
for all £>3. 

Proof. We will first characterize the vector space K[x ly ... , x n ] j {p n ,k) K ^ ■ Observe that the coefficient 
matrix A := (v ; ) ;e [ n _ fc ] is actually a Vandermonde matrix and in particular every square submatrix 
of A is invertible. Furthermore, the polynomials v ; - with j e [n — k] are homogeneous of degree 
one. Thus, for any k variables of ...,x n }, without loss of generality say, x-y, x k , we have that 
{xi, x k , Vi, . . . , v n _ k } is a basis for the homogeneous polynomials of degree one. This is just another 
way of saying that removing the columns belonging to x x , . . . , x k from A, the resulting square submatrix 
is invertible. It follows 



K[x 1 ,...,x n ]l (Fn,k) K[ x] ~ K[x 1 ,...,x k ,v 1 ,...,v n _ k ]j (F n ,k} K [x] 

= K[x 1 ,...,x k ]l(Ti 3 ) Km . 



As the substitution preserves degrees, homogeneity, etc., it follows, that any degree-compatible order 
ideal has to have \0 =1 \ = k, G =2 = {xy | x,y e G= 1 }, and =l = for all I > 3. 

The other direction follows immediately as each order ideal with \0 =x \ = k, =2 = {xy i,y e 
=1 }, and =l = for all I > 3 is actually a degree-compatible order ideal such that K[xi, ■ ■ ■ , x n ] = 
{ F n,k) K [ X j ® {&)k as vector spaces, by the argumentation above. □ 

Note that the order ideals of F n k indeed correspond to the fc-cliques of the complete graph on n 
vertices: If e A(F n k ), then =1 = {x ti , . . . ,x ik } and x^x^ e =2 if and only if X;.,^ e =1 . If we 
now remove all elements of the form x 2 with x,- e =1 , and there are k of those, then 

lj 1; 

fc(fc-l) 

\0= 2 \{x 2 \ Xij Z0=i}\ = ^^, 



the size of a fc-clique. We are ready to state the main result of this section Consider the following 
problem: 

Max-Order ideal of ideal: Let M c K[X] be a system of polynomials generating a zero-dimensional 
ideal and let c e Z T " be a preference on the monomials. Compute an order ideal supporting an 
^■-border basis of (M) K [ X ] with maximum score with respect to c. 
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Theorem 5.3. IMax-Order ideal of idealI ls NP-hard. 



Proof. The proof is by a reduction from the NP-hard ?c-Clique along the lines of the proof of Theo- 
rem 5.1 Let r = (V, E) be a graph with n := |V| and k e [n] be an instance of ?c-Clique We consider 
M := F n k and define c e Z T <3 via 

fl, if m = x u x v and (u, v) G E; 
'" 1 0, otherwise. 



By Lemma 5.2 we have that the degree-compatible order ideals of {M) K ^ are in one-to-one correspon- 
dence with the fc-cliques of the complete graph on n vertices. Similarly to the proof of Theorem |5.1[ the 
score of an order ideal is just the number of edges between the corresponding k vertices of the graph, 
so the maximum score is k(k — l)/2 if and only if the graph contains a fc-clique. Thus, we obtain a test 
whether r contains a clique of size k. □ 

6. Computational results 

We performed a few computational tests to verify the practical feasibility of our method. All compu- 
tations were performed with CoCoA 4.7.5 (H7J) and scip 1.1.0 ((ED) on a 2 GHz Dual Core Intel 
machine with 2 GB of main memory. We performed computations on various sets of systems of polyno- 
mial equations. The employed methodology was as follows. We first computed a border basis using the 
classical border basis algorithm. From the last run of the algorithm we extracted the L -stabilized span 
and brought it into canonical form. We generated the constraints fl3.7] ) from the order ideal that we 
obtained; from the L -stabilized span in matrix from, we generated the constraints fl3.8ty . As we needed 
to get access to the L -stable span computed in the last round of the border basis algorithm, we had 
to use an implementation of the border basis algorithm in CoCoA-L which is slower in terms of speed 
compared to a C or C++ implementation. We then transcripted these constraints into the CPLEX LP 
format which served as input for scip. For the optimization we chose various preference functions. 
One was a random function, and the other one was chosen with the intent to make the optimization 
particularly hard by giving monomials deep in the order ideal negative weights and assigning positive 
weights for the outer elements. 

In all cases the optimization, i.e., the computation of the weight optimal order ideal was performed 
in less than a second, whereas the actual calculation of the initial border bases was significantly more 
time consuming. As indicated before, this is not unexpected as the computation of the L -stable span 
is significantly more involved than computing a weight optimal order ideal (in the worst case double 
exponential vs. single exponential). When computationally feasible we also counted all feasible order 
ideals with scip, which basically means enumerating all feasible solutions. This in fact is equivalent 
to optimizing all potential preference functions simultaneously and thus emphasizes the computational 
feasibility. 

We considered 7 systems of polynomials of various complexity in terms of number of variables and 
order ideal degree. We would have liked to test significantly larger instances but we were not able to 
compute the initial border basis (or more precisely the L-stable span), neither with our implementation 
nor with the C/C++ implementation of CoCoA 5 (BBasis5). This shows once more that the limiting 
factor is the actual computation of the L-stable span. 

In Table [I] we report our results. 

7. Concluding remarks 

We provided a way to characterize all degree-compatible order ideals that support a border basis 
for a given zero-dimensional ideal / by borrowing from combinatorial optimization and in particular 
polyhedral theory. We established a one-to-one correspondence of the integral points of a certain 
polytope, the order ideal polytope, and those degree-compatible order ideals that support a border 
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polynomial system 


order ideal signature 


optimization [s] 


counting [s] 


# order ideals 


■5 2 3 

x ,xy +y 


(1,3,1,1,1) 


< 0.01 


0.U2 


3 


vanishing ideal of the 

points (0,0, 0,1), 
C1,0,0,2J, (3,0, 0,23, 
(5,0, 0,33, (-1,0,0,43, 
(4,4,4,5), (0,0,7,6)). 


(1,4,2) 


< U.Ul 


(J.U2 


45 


x + y + z — u — v, 


(1,4,5) 


< 0.01 


0.35 


1,260 


222 
x —x,y -y,z -z, 

2 2 

u — U, V — V 










x + y+z — u — v, 


(1,4, 7,6) 


0.02 


51.50 


106,820 


3 3 2 

x —x,y -y,z -z, 

2 2 
u — u, v — V 










x + y+ z — u — v, 


(1,4,8,9) 


0.U2 


53. uu 


108, 9UU 


333 
x -x,y -y,z -z, 

2 2 
u — u, v — V 










x + y+ z — u — v, 


(1,4,9,12, 9) 


U.Uo 


3(J(J.(J(J 


> 1,349,154 


333 
x J -x, y -y, -z, 

3 2 
ir - u, v z - v 










x + y + z — u — v + a, 


(1,5,9) 


< 0.01 


8.68 


30,030 


222 
x z -x, y -y, z z -z, 










222 
u — u, v — v, a — a 











Table 1. Computational results. The first column contains the considered polyno- 
mial system. The second column contains the degree vector of the order ideal, i.e., 
(dim I- 1 jl- 1-1 ). starting with i = 1 and J-° := 0. The third column contains the av- 
erage time (in seconds) needed to optimize a random preference over the order ideal 
polytope (we performed 20 runs for each system). The fourth column contains the time 
(in seconds) needed to count all admissible degree-compatible order ideals and the last 
column contains the actual number of admissible degree-compatible order ideals. The 
'*' indicates that the counting had been stopped after 300 seconds. The number of order 
ideals reported in this case is the number that have been counted up to that point in 
time. 



basis of I. This connection in particular links the ideals to their combinatorial structure of the factor 
spaces. 

Using this polytope we adapted the classical border basis algorithm in order to be able to compute 
border basis for general degree-compatible order ideal based on a preference ordering on terms con- 
tained therein. Effectively, the algorithm can be used for any integral point contained in the order ideal 
polytope. 

We also showed that computing a border basis for a preference on the monomials one might want to 
have included in the order ideal is NP-hard and thus we cannot expect to be able to efficiently compute 
preferred order ideals in general although it is merely a basis transformation. On the other hand, 
this is not restricting the applicability of our method in any practical application because the preceding 
computation of the L -stable span dominates in terms of computational complexity. We finally presented 
a few computational results showing the applicability of our method for actual computations. 
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